A class of exactly solvable pairing models 
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We present three classes of exactly solvable models for fermion and boson systems, based on the 
pairing interaction. These models are solvable in any dimension. As an example we show the first 
results for fermion interacting with repulsive pairing forces in a two dimensional square lattice. 
Inspite of the repulsive pairing force the exact results show attractive pair correlations. 
PACS number: 71.10.Li, 74.20.Fg 



en ■ 

<n: 

o- 

o ; 
u • 

Of 

+-* . 



-a 
c 
o 
o 



> 

o 



c3 



i 

C 

o 
o 



X 



Exactly solvable models have played an important role 
in understanding the physics of the quantum many body 
problem, especially in cases where the system is strongly 
correlated. Such situations arises e.g. in one dimensional 
(ID) systems of interest for condensed matter physics 
and also in strongly correlated finite fermion systems as 
atomic nuclei. In both branches of physics the study of 
exactly solvable models has been pursued since long with 
enormous success. 

In ID quantum physics, the exactly solvable models 
can be classified into three families. The first family be- 
gun with Bcthe's exact solution of the Heisenberg model. 
Since then a wide variety ID models has been solved us- 
ing the Bethe ansatz (for a recent review see H). A 
second family of models are the so called Tomonaga- 
Luttinger models [|l| which are solved by bosonization 
techniques and which revealed non Fermi liquid proper- 
ties in ID. These systems are now called Luttinger liq- 
uids. The third family, proposed by Calogero and Suther- 
land, are models with long range interactions. They have 
been applied to several problems 0] like spin systems, the 
quantum Hall effect, random matrix theory, etc... 

Several exactly solvable models have been developed 
in the field of nuclear physics from a different perspec- 
tive 0]. In these models the hamiltonian is written as a 
linear combination of the Casimir operators of a group 
decomposition chain ideally representing the properties 
of a particular nuclear phase. Typical examples are the 
Elliot SU(3) model describing nuclear deformations and 
rotations and the U(6) Interacting Boson Model with 
its three dynamical symmetry limits describing rotational 
nuclei (SU(3)), vibrational nuclei (U(5)) and gamma un- 
stable nuclei (0(6)). These models were extremely useful 
in providing a simple understanding of some prototypical 
nuclei. 

The impact of the exactly solvable models in condensed 
matter physics and in nuclear physics is so enormous that 
one hardly can believe that the exact solution of the Pair- 
ing Model (PM), of great interest for both fields, passed 
almost unnoticed till very recently Q. The PM, was 
solved exactly by Richardson in a series of papers in the 
sixties §-§. 

Independently of Richardson's exact solution, it was 
recently demonstrated || that the PM is an integrable 
model. The pairing model may turn out particularly 



interesting, since recent work H has shown that the 
pure repulsive pairing Hamiltonian in a 2D lattice can 
be solved exactly in the thermodynamic limit revealing 
strong superconducting fluctuations. The importance of 
this finding stems, of course, from the fact that high T c 
superconductors apparently acquire their superconduct- 
ing properties through the repulsive Coulomb interac- 
tion. 

We will derive in this letter three families of exactly 
solvable models based on the pairing interaction for 
fermion systems as well as for boson systems. The most 
important feature of the new set models we will present 
is that they are exactly solvable in any dimension. In 
pO| we have advanced a numerical solution for a three 
dimensional confined boson systems, here we will give 
preliminary results for a fermion system in a 2D lattice. 

Since the proof of integrability and the derivation of 
the exact solutions is completely analogous for fermions 
and bosons, we will develope both systems in parallel. 
In what follows whenever there are different signs, the 
upper one will correspond to bosons while the lower one 
to fermions, and we will refer indistinctly to bosons and 
fermions as particles. 

Let us begin our derivation by defining the three oper- 
ators 



ni 



5>L«n» , 4=E°L«k = (^) t a) 



which close the commutator algebra 



n, , 



= 2S W A\ , \A h A\, 



= 2S w (Qi±2rii) (2) 



In eq. (H) the pair operator Aj creates a pair of parti- 
cles in time reversal states with a' (a) the particle creation 
(annihilation) operator and f2; being the degeneracy of 
level I. 

The number operator n; and the pair operators A\ , Aj 
in each level I close the commutator algebra of the groups 
SU(2) for fermions or SU(1, 1) for boson systems 

The three generators of these algebras can be written 
in terms of the previously defined pair and number oper- 
ators as Kf = \ni ± ifi;, and K+ = \A\ = (iff)*. 

These generators obey the more familiar commutation 
relations of the SU(1, 1) and SU{2) group algebras 



[Kl K+] = 5 W K+ , [K+, Kj] = T26 IV K? (3) 

The difference between the algebra of SU(1, 1) and SU(2) 
appears in the sign of the second commutator. 

The Hilbert space of N particles moving in L single 
particle levels can be classified according to the product 

■ ■ x SU (2) L for fermions 
x SU (1, l) L for bosons. 
A complete set of states in this Hilbert, which are space 
spanned by the pair operators (pi), is given by 



of groups SU (2) 1 x SU (2) 2 
or SU(1, 1) 1 x SU (1,1) 2 x 



\ni,n 2 ,---,n L ,i/) 
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j^m j^n 2 



■4 ni 



(4) 



where TV is a normalization constant. The possible num- 
ber of pairs in each level is < ni < J7;/2 for fermion 
systems or < ni < N/2 for boson systems. A state 
\v) = \v\V2 • • • vl) of unpaired particles is defined as 



Ai \v) = , ni \u) = i>i \v) 



K?\» 



where d, = (%■ ± M 1 ) and A = 2M 



\")=di\v) (5) 

M being the 

number of pairs and v the total number of unpaired par- 
ticles. We will borrow from Nuclear Physics the name 
Seniority for the number of unpaired particles in each 
level v\. 

In the product spaces mentioned above a model is in- 
tegrable if there are L independent global operators com- 
muting with one another. These operators are the quan- 
tum invariants and their eigenvalues are the constants of 
motion of the system. In looking for these operators, let 
us define the most general combination of one and two 
body operators in terms of the K generators with the 
condition of being hermitian and number conserving: 



Ri = K? 




K K v) 



ty w k?k?,} 



(6) 



Up to now the matrices X and Y are completely free, 
but we will fix them imposing the condition that the R 
operators should commute among one another to define 
an integrable model. The condition [Ri,Rir] = will be 
fulfilled if they are antisymmetric and satisfy the follow- 
ing equation 



Y t jX 



jk 



Y kt X 



ki-A-jk 



XkiXij — 



(7) 



An analogous condition has been encountered by 
Gaudin |11| in a spin model known at present as the 
Gaudin magnet. His model is based on R operators simi- 
lar to (||) but without the one body term. He found three 
different solutions for the condition (M) which can be 
grouped together in compact form as Xw = s ., ( 7 _ ,, , 
Yw = 7 cot [7 (r]i — ni/)}, where the different models are 
distinguished by the value of 7: 

I. The rational model: 7 — > 



Xn, = Y w = 



Vi - Vi' 
II. The trigonometric model: 7=1 



X h 



1 



, Y w = cot (iji -rji') 



sin (ni - rji' ) 
III. The hyperbolic model: 7 = 1 

1 



(8) 



(9) 



Xw = 



sinh(?7; - ni>) 



, Yw =coth (rji-r)i>) (10) 



rji is an arbitrary set of non-equal real numbers. Each 
solution gives rise to an integrable model and any combi- 
nation of the R operators produces an integrable hamil- 
tonian. Since the three models have quite different prop- 
erties, we prefer to continue our derivation in an indepen- 
dent way. It is worthwile to mention here that if we relax 
the condition of number conservation on the set of oper- 
ators (pf) there is a general solution in terms of elliptic 
functions |TjJ . 

The rational model has been recently proposed in ref. 
H to demonstrate the integrability of the PM hamilto- 
nian. Indeed the PM hamiltonian can be obtained from 
the rational model by means of linear combination of R 
operators Hp = 2^ ; ViRJ plus an appropriate constant 
to give 



tfp = E^ + ?E^ 



(ii) 



where the free parameters rji have been replaced by the 
single particle energies ei . 

Then, the PM hamiltonian is diagonal in the basis 
of common eigenstates of the R operators. But, also, 
any function of the R operators defines a valid integrable 
hamiltonian. If we want to restrict ourselves to hamil- 
tonians containing at most one and two body terms, the 
most general linear combination of R operators is 

H = 2 J2 ZiKi = 2 E e ' K i + 2 9 Yl ( £l ~ ev ) Xw K t K v 

T2gY,(zi-ei') Y u' K ? K ° (12) 

In (IT2]) the X and Y matrices can be given formally by 
one of the three models (PHlO|), and for each of the models 
they are functions of an arbitrary set of r/'s and a pairing 
strength g. Moreover the set of e's in the linear combi- 
nation (|l2j) is also arbitrary, so that the total number of 
free parameters in defining an integrable hamiltonian in 
each of the three models is 2L + 1 for fermion as well as 
for boson systems. 

We would like to emphasize that the demonstration 
of integrability given above for the three models does 
not imply that they are exactly solvable, i.e. that the 
complete set of common eigenstates of the R operators 
can be obtained. As mentioned before, Richardson only 



worked out the eigenstates for the PM hamiltonian and, 
as a matter of fact, the authors of ref. || were not aware 
of the existence of that exact solution. Very recently |1J] 
the eigenvalues of the R operators of the rational model 
for fermion systems were obtained, within a conformal 
field theory formalism, in the fully paired subspace (V/ = 
0) and with constant degeneracy Hi — 2 . 

To begin our derivation we will propose an ansatz for 
the eigenstates of the R operators in the Hilbert space of 
states ([|) which is a generalization of the ansatz used by 
Richardson to find the eigenstates of the PM. The exact 
eigenstates for the three models can be written as prod- 
uct pair wavefunction acting on the space of unpaired 
particles \v) 



m 






M 

\[Bl\v), Bl = ^ Ul (e a ) Kt 1181 

a=\ I 



The function u depending on the pair energies e, has 
the form of the eigenstate of the one pair problem in each 
of the three models and the set of parameters e, are left 
as free parameters to be fixed in order to fulfil the L 
eigenvalue equations 



Ri |*> = r< Itt) 



(14) 



The collective amplitudes in the pair operators B for 
each model are 



uii = 1H (e Q ) 



1 



2th 



ii-iii 



Ui (e«) 



1 



sn (e a - rji) 



(15) 



(16) 



We will treat on equal footing the trigonometric and 
the hyperbolic models to solve the eigenvalue equation 
( |l4| ) . To embody both derivations in the same formalism 
we will use the symbols sn for sin or sinh, cs for cos 
or cosh and ct for cot or coth, not to be confused with 
elliptic functions. Note that we have already used sn in 

©• 

Here we summarize the final results for the three mod- 
els leaving the details of the derivation for a future pub- 
lication. 

Model I 
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Models II and III 



-^ T 4gV X 
f]i - Vj ^Vi - e c 



(17) 



(18) 



1 =F 2 g£ dj ct(e a - rjj) ± 2g £ ct(ep — e a ) = (19) 




FIG. 1. Real and positive imaginary parts of the pair en- 
ergies e a for a 6x6 lattice at half filling as a function of g. 



r i = d l <lT 2g 






(20) 



In order to obtain the pair energies e a , given a set of 
parameters r\ and a pairing strength g, one has to solve 
the coupled set of M nonlinear equations ([17]) for the ra- 
tional model or ( |l9| ) for the trigonometric or hyperbolic 
models respectively. In the limit g — > one imediately re- 
alises that the equations ([17], |l^) can only be fulfileed for 
e a — > 2?7i. Then the amplitudes u a i in ( J15] , Hq) become 
diagonal and we see that the states ( |13|) are equivalent 
in this limit to the complete set of uncorrelated states 
(j|). From this g — + limit one can construct the ground 
state (with pairs filling the lowest possible states), and 
the configurations of the succesive excited states. For 
example, the first excited state is obtained by promoting 
the highest energy pair to the next upper 2g value, or 
by breaking a pair (removing a pair energy) into two un- 
paired particles (increasing the senioriy v by two). One 
then follows the trajectory of the pair energies e a for 
each of the configurations as a function of g solving the 
equations ( |17] |19|) . 

For bosons systems the pair energies stay always real, 
but for fermion systems the pair energies can be either 
real or complex conjugate pairs. In the latter case there 
might arise singularities in the solution of the equations 
(O, h9) for some critical value of the pairing strength 
g c for which two or more pair energies acquire the same 
value. It can be shown pi that each one of these critical 
g values is related to a single particle level i and that 



TABLE I. 

6x6 lattice. 


Single particle energies and degeneracies for the 


*i -4 


-3-2-1 12 3 4 


Qj 2 


8 8 8 20 8882 



TABLE II. Correlation energy densities and Interaction 
energy densities for various values of g for the 6x6 lattice 
at half filling. 



at the critical point 1 — 2di pair energies must be equal 
to 2r]i. These sigularities cancel out in the calculations 
of energies which do not show any discontimuity in the 
vicinity of the critical points. 

The eigenvalues of the R operators, given by (|l^) or 
( P0| ) respectively, are always real since the pair energies 
are real or complex conjugate pairs. Each solution of the 
nonlinear set of equations produces an eigenstate com- 
mon to all Ri operators, and consequently to any hamil- 
tonian written as linear combination of them. The cor- 
responding eigenvalue is the linear combination of the 
r, eigenvalues. As mentioned before, the most impor- 
tant feature of the three families of models is that they 
are exactly solvable in any dimension. The dimension- 
ality enters through the degeneracies di = vi/2 — fij/4. 
Also through the same term enters the information about 
states with nonzero seniority (broken pairs). Symme- 
try breaking terms, like anisotropic hoppings or disor- 
der, might lift up the degeneracies but still leaving the 
hamiltonian as exactly solvable. In such cases the dimen- 
sionality would show up in the density of states. 

Next we will present the first results for the ratio- 
nal model of fermions in a 2D square lattice of P x P 
sites with periodic boundary conditions and a repul- 
sive pairing interaction. Assuming a restricted hopping 
term between nearest neighbors, the single fermion lev- 
els are £& = — 2 (cos k x + cos k y ), with k a = 2Trn a /P and 
—P/2 < n a < P/2. Here we will consider a 6 x 6 lat- 
tice at half filling (M — 18) with a PM hamiltonian for 
which r\i — Si in (|17|), but the properties we will discuss 
are quite general and independent of the latter choice. 

The single fermion energies £& and the corresponding 
degeneracies $\- are displayed in the Table I. In the limit 
g = the groundstate is obtained by distributing the 
M = 18 pairs in the lowest possible states. 

In Fig. 1 we show the real and imaginary part of the 
pair energies. Only real part and the positive imaginary 
part are shown for each complex conjugate pair. The 
first level accomodates one pair which is forced to be 
real. The next three levels accomodate four pairs each, 
forming two pairs of complex conjugates in each level. In 
the last level we can put five pairs of particles forming two 
pairs of complex conjugates and one real pair energy. For 
the critical value g\ = 0.1708 the first five pairs become 
equal to 2e2 = — 6 as discussed before. As seen in Fig. 
1, the real part of the two first complex conjugated pairs 
together with the first real pair energies cross at —6 while 
the imaginary parts go to zero. A similar situation arises 
for the second critical point at gi = 0.2407 in which five 
pairs become equal to 2e 3 = —4. 

In spite of the quite involved behavior of the pair en- 
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ergies, the total energy displays a smooth behavior as a 
function of g. In Table II we give some values of the cor- 
relation energy density ((E(g) — EHF(g))/P 2 ) and the 
interaction energy density ((E(g) — E(Q))/P 2 ) as a func- 
tion of g. The interaction energy increases with increas- 
ing values of g, showing a tendency to saturation that is 
consistent with its vanishing in the thermodynamic limit 
|l4f . The correlation energy density decreases almost lin- 
early with g. These attractive pairing correlations based 
on a repulsive interaction give further numerical support 
to the work of Shastry 19] who found quasi-long ranged 
order in the thermodynamic limit of the model at half 
filling. 

In summary, we have presented three families of new 
exactly solvable models based on the pairing interaction 
for fermion and boson systems. These models have the 
important feature of being solvable in any dimension. We 
have presented preliminary results for the properties of 
the exact solution of the rational model in a 2D square 
lattice with repulsive pairing interaction. This model 
may be useful to study features of high Tc superconduc- 
tivity because, inspite of the purely repulsive bare force, 
the exact solution shows attractive pair correlations. 

After completing this work we learned of a recent 
preprint |15| in which the hyperbolic model is presented 
for electrons in ultrasmall superconducting grains, how- 
ever without applications. This model is equivalent to 
our models II and III for fermions in ID with seniority 
zero. 
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